The Effects of Finite Emission Height in Precision Pulsar Timing 



Caroline D'AngekQ 
Department of Astronomy and Astrophysics, University of Toronto, 
50 St. George Street, Toronto, ON M5S 3H4, Canad£\ 

Roman R. Rafiko\[l 

CITA, McLennan Physics Labs, 60 St. George St., 
University of Toronto, Toronto, ON M5S 3H8, Canada 
(Dated: February 5, 2008) 

Precision timing is the key ingredient of ongoing pulsar-based gravitational wave searches and tests 
of general relativity using binary pulsars. The conventional approach to timing explicitly assumes 
that the radio emitting region is located at the center of the pulsar, while polarimetric observations 
suggest that radio emission is in fact produced at altitudes ranging from tens to thousands of 
kilometers above the neutron star surface. Here we present a calculation of the effects of finite 
emission height on the timing of binary pulsars using a simple model for the emitting region geometry. 
Finite height of emission h changes the propagation path of radio photons through the binary and 
gives rise to a large spin velocity of the emission region co-rotating with the neutron star. Under 
favorable conditions these two effects introduce corrections to the conventional time delays at the 
microsecond level (for a millisecond pulsar in a double neutron star binary with a period of several 
hours and assuming h — 100 km). Exploiting the dependence of the emission height on frequency 
(radius-to-frequency mapping) and using mult i- frequency observations one should be able to detect 
these timing corrections even though they are formally degenerate with conventional time delays. 
Although even in the most accurately timed systems the magnitude of the finite emission height 
effects is currently somewhat below timing precision, longer-term observations and future facilities 
like SKA will make measurement of these effects possible, providing an independent check of existing 
emission height estimates. 

PACS numbers: 95.30.Sf - 04.25.Nx - 04.80.-y 



I. INTRODUCTION 

Since their discovery radio pulsars have been recog- 
nized as precise time-keepers. This unique property has 
been extensively exploited to test general relativity since 
1975, when Hulse and Taylor discovered the first dou- 
ble neutron star (DNS) system PSR 1913+16 0. Thirty 
years later, pulsar surveys have found seven other DNS 
systems (among them a double pulsar system J0737- 
3039, in which both neutron stars are seen as radio pul- 
sars, see 0) and many white dwarf- neutron star (WD- 
NS) binaries. These systems have been used not only 
to test a number of general relativistic effects but also 
to set limits on the stochastic gravitational wave back- 
ground 0- 

Precise monitoring of the times of arrival (TOAs) of 
pulses on Earth is the key to success in this game. 
Although timing accuracy at the several /zs level has 
been routinely achieved in a number of systems, only 
three WD-NS binary pulsars - J0737-4715, J1713+0747, 
J1909-3744 HH " have been timed to 0.1 - 0.3 ^s preci- 
sion. At the same time, current instruments have not yet 
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reached the level of intrinsic uncertainty present in pul- 
sars (which will ultimately set a limit on their precision, 
see (5()- Ambitious new projects like the Parkes Pul- 
sar Timing Array (PPTA), which aims to achieve daily 
timing residuals < 0.1 /is for more than ten millisecond 
pulsars @, and proposed facilities like the Square Kilo- 
meter Array [SKA; 5f, demand an understanding of the 
systematic timing uncertainties at the level of 1 ns 0. 
With these new instruments, effects previously hidden in 
error will become observable and some of them may have 
not been accounted for by the current timing models. 

The conventional physical model for interpreting tim- 
ing data in terms of the binary system's parameters, 
which was advanced by Damour and Deruelle la, here- 
after DD] and extended by Damour and Taylor |9j, here- 
after DT] , treats the binary as two point masses in orbit 
around each other and assumes the radio emission re- 
gion to be located at the very center of the radio pulsar. 
However, theoretical models of pulsar radio emission, de- 
tailed analysis of pulse structure and polarimetry of iso- 
lated pulsars suggest that their radio emission is in fact 
produced high in the magnetosphere. In slow pulsars it 
occurs at typical distances of order several hundred (or 
thousand in some cases) kilometers from the neutron star 
[ToL [Til n"3, IT^ . For millisecond pulsars, which are more 
appropriate for precision timing, the typical height of the 
emission re gion h is estimated to be several tens of kilo- 
meters ^10, although at low frequencies h may be as 
high as « 100 km. 
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A finite emission height has twofold effect on pulsar 
timing. First, the displacement of the emission region 
from the neutron star center affects photon propagation 
through the curved space-time within the binary. Second, 
a high-altitude emission region of a rapidly spinning mil- 
lisecond pulsar must be moving at a rather large speed 
(in some cases reaching 0.2c [lq . where c is the speed 
of light), leading to severe aberration of the pulsar radio 
beam. For a given spin period of a pulsar P p — 2ir/Q p 
the upper limit on the emission height h is set by the 
size of the light cylinder r LC = c/Cl p » 500(P p /10 ms) 
km, while the lower limit on h is set by the neutron star 
radius Rns ~ 10 km. The latter also sets a lower limit of 
QpRNS sin a 0.02c(P p /10 ms) -1 sina on the rotational 
velocity of the emitting region, where a is the angle be- 
tween the magnetic and spin axes of the pulsar. 

Given the rapidly improving accuracy of pulsar timing 
and the increasing demand to understand systematic tim- 
ing effects at the nanosecond level, which are driven by 
ongoing and future projects, the aforementioned effects of 
a non-zero emission height merit closer inspection. The 
aim of this paper is to extend the DT model of pulsar 
timing to the case where the emission comes from a re- 
gion high in the magnetosphere, and determine whether 
this physical change produces any measurable results. 

An outline of the paper runs as follows. In sec. ITT1 we 
describe the geometry of the emission in a dipolar mag- 
netic field. In sec. IIIII we consider different timing effects 
and compute corrections that result from considering a 
finite emission height. We discuss the measurability of 
these corrections in sec. IIVI and limitations and possible 
refinements of the model in sec. [V] 



II. GEOMETRICAL MODEL OF RADIO 
EMISSION 



We consider a binary system with a semimajor axis 
a, eccentricity e, and inclination i with respect to the 
observer's line of sight. The instantaneous position of the 
pulsar is characterized by the true anomaly ip measured 
from the ascending node of the pulsar orbit, and u is the 
longitude of periastron, see Figure We use R p and R c 
for the position vectors of the pulsar and its companion 
with respect to the barycenter of the binary. In making 
numerical estimates we will adopt a fiducial model of the 
pulsar binary (motivated by the properties of the double 
pulsar J0737-3039, see 0) which consists of two neutron 
stars (M p = M c = 1.4 Mq are the assumed pulsar and 
companion masses) in a compact orbit with a semimajor 
axis a — 10 11 cm (orbital period Pi, = 2.9 hrs). We take 
pulsar spin period to be P p = 10 ms and assume the 
height of the emission region to be /17 = 10 7 cm. 

To illustrate our calculations we adopt a specific model 
for the geometry of pulsar radio emission closely related 
to the conventional rotating vector model [RVM; 0] . In 
the comoving coordinate system of the pulsar we assume 
neutron star magnetic field to have a dipolar geometry 
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FIG. 1: Orbital configuration and spin orientation of the pul- 
sar. Here s p is the spin axis for the neutron star, no is the 
direction to the observer, rh is the magnetic field axis, R p 
is the position vector of the pulsar. We use the orientation 
notation of DT so that the orbital angular momentum vector 
L constitutes angle i with —no. See text for details. 



(which should be a good approximation at the distances 
where radio emission is produced), with the magnetic 
axis rh (|rh| = 1) making an angle a with the pulsar 
spin axis s p (|s p | = 1), see Figure |2] Vector s p makes an 
angle C with the unit vector no from the binary system to 
observer at Earth (s p ■ no = cos£), while the projection 
of s p on the plane of the sky makes an angle 77 with the 
ascending node of the binary orbit, see Figure ^ While 
the vector s p is fixed, the orientation of rh changes in 
time because the pulsar rotation always maintains the 
relation 



cos a. 



(1) 



The position of the emitting region with respect to the 
pulsar center is hr, where |f| = 1 and h is the emis- 
sion height, see Figure f3J Emission is assumed to oc- 
cur strictly along the local B field since it is likely pro- 
duced by the highly relativistic plasma streaming along 
the open field lines, so that the direction of emission d 
(|d| = 1) is related to f and m as 



d = 



3 cos 9r — rh 



(3cos 2 0+l)V2' 
where 9 is the angle between f and rh, i.e. 



m • r = cos ( 



(2) 



(3) 



Vector d makes an angle p with the magnetic axis rh (see 
Figure fSJ such that [see eq. ©] 



cosp 



3 cos 2 6-1 
(3cos 2 + l)V2- 



(4) 
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FIG. 2: Model of the emission cone geometry of the pulsar. Here we show the direction of emission d when d = no and the 
pulse is observed. Note that at the moment of emission rh has two values, which correspond to the two crossings of the emission 
cone by the observer's line of sight. The coordinate system (x, y, z) is also displayed. 



We will be considering emission by all sources at the 
same height h belonging to a single flux surface for which 
the condition ® with a given value of 9 is fulfilled. Thus, 
h and 8 are fixed for this calculation (but see sec. IV Al 
for the consequences of relaxing this assumption). All 
such sources form a circular emission region around vec- 
tor rh, which gives rise to an emission cone with an open- 
ing angle 2p rigidly connected to (and co-axial with) the 
magnetic axis. When the emission cone crosses our line 
of sight as a result of pulsar rotation (that is, for some 
point hr in the emission region, no becomes parallel to 
d) a pulse of emission is sent to the observer on Earth. 
There are two emission episodes per spin period of pulsar 
corresponding to the two emission cone crossings by no, 
separated by 2<E>o in pulse phase (see Figure |2J, where 



cos $o 



cos p — cos C cos a 
sin £ sin a 



(5) 



By measuring the TOAs of these pulses one can obtain 
information about the motion of the pulsar and the prop- 
agation effects inside the binary. 

The emission cone crosses our line of sight and makes 
the pulsar observable only if |£ — a\ < p. In the absence 
of aberration caused by the motion of the emission region 
(i.e. when pulsar is not spinning and is at rest with re- 
spect to observer on Earth), favorable conditions for the 
detection of this conal emission, 



d = n , 



(6) 



can be achieved only for some specific orientations of vec- 
tors f and rh, which we call fo and rho. To find these 
zeroth order (or static) values of f and rh we form a sta- 
tionary orthogonal basis using s p and fio- 



(x,y,z) 



s„ x n s p — cos (ho 



x n Q 



x n | 



-,n 



(7) 



FIG. 3: An illustration of the dipole model used in our analysis. The field lines emanate from the surface of the star about the 
magnetic field axis, rh, and the emission is produced at the point hf , in the direction d, which is tangent to the field line. 



Then, using equations , J2J , 10 , and © one can easily 
find that 

cos a — cos C cos p 

mo = ± sin $o sin a x H y (8) 

sin( 

+ cos p z, 

sin $o sin a _ cos a — cos C cos p „ 

r ° = ± ^ a — x + 5 a~- — 7 y ( 9 

3 cos 6 3 cos 6 sin Q 

cos p + [3 cos 2 + 1] 2 ^ 
3cos# 

where sin $o can be found from equation JSJ . The ± 
ambiguity in the x-components of mo and f o corresponds 
to the two crossings of the emission cone by no which 
results in two possible orientations of rh and f , see 
Figure |21 The upper and lower signs are for the leading 
and trailing pulses, respectively (cone crossings). 

Slow pulsars typically have narrow beams meaning 
that p < 1, $ < 1, and 9 w (2/3) 1/2 p < 1. In 
this case it follows from equations 10 and I|1U|I that 
mo, X) &o, y , fo,x, ?o,y ~S< 1 (remember that |C — «| < p) 
so that both mo and fo are almost aligned with no. 
On the other hand, millisecond pulsars (the objects 
most suitable for precise timing), have wide beams with 
p, $o, @ ~ 1> so that their mo and fo are generally quite 
mis-aligned with no. 

When the emission region is moving, the concomitant 
aberration changes the emission direction d as seen by 
an observer at rest, so that the radio pulses are detected 



only at f and rh, which are somewhat different from fo 
and riio. The conventional timing model (DD) assumes 
that the motion of the emission region comes only from 
the orbital motion of pulsar with speed v D = (3 c (the 
proper motion of the binary as a whole is undetectable 
through timing). A finite height of emission introduces 
the spin velocity v s = (3 s c, since the emission region co- 
rotates with the pulsar. In the case of millisecond pulsars 
this spin velocity is a considerable fraction of the speed 
of light ^(| even if emission takes place at the surface 
of the neutron star (v s w 0.2c for an orthogonal rotator 
with spin period 1 ms and h ~ 10 km, which is much 
larger than v Q ~ 10~ 3 c). 

In the general case of nonzero orbital and spin veloci- 
ties of the emitting region unit vectors f and rh can be 
written as 

f = f + A r + A s r + A so r, (10) 
rh = rh + A G m + A s m + A so m, (11) 

where A Q denotes corrections due to v Q (which have 
been previously calculated in the framework of DD tim- 
ing model), A s corresponds to the corrections due to 
v s only, while A so corresponds to the coupling be- 
tween v c and v s . The derivation of all these cor- 
rections is presented in Appendix and the expressions 
for A D r, A s r, A so r, A D m, A s m, A so m are given by equa- 
tions 1|A4|) - (|A5|) and HA12MA14[I . For tne parameters of 
our fiducial pulsar binary model the typical relative mag- 
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nitudes of the different correction terms listed in equa- 
tions (TDJ and (fTTf) are 



orbit. This delay can be expressed as: 



V ° ~ m-3 -1/2 

£o = — w 10 a n , 
c 

£s = ^ = 0.21 fc 7 (JJL. 

c V 10 ms 



(12) 
(13) 



e so = e s s » 2 x 10- 4 /i 7 a n 1/2 f t^M , (14) 

\ 10 ms / 

correspondingly, where hi = /i/(10 7 cm) and an = 
a/(10 n cm). 



III. EFFECT ON TIMING 

We are now in a position to calculate the corrections 
to the time delays that occur in binary pulsar systems. 
These delays appear when one relates the observer's 
proper TOA of a given radio pulse r to the proper time 
of its emission in the comoving frame of the pulsar T e . As 
demonstrated by DD this transformation can be done in 
several steps, going first from r a to the coordinate time 
of arrival t a (measured in a frame comoving with the bi- 
nary barycenter), which is linked to the coordinate time 
of emission t e through the two propagation delays - the 
Romer (Ar) and Shapiro (Ag) delays (see eq. [4] of DD). 
In its turn, the coordinate time t e is directly related to 
the proper emission time in the pulsar frame T e via the 
Einstein delay A^, so that t e = T e + Ae- Finally, the de- 
viation of the proper time of emission T e from the purely 
periodic pattern in the pulsar frame is determined by the 
so-called aberration delay A^. 

The introduction of the nonzero emission height modi- 
fies some of these delays. The Romer and Shapiro delays 
get affected because of the change in the path of the radio 
beam. The aberration delay gets modified because finite 
height introduces spin motion of the emission region (in 
addition to the orbital motion) which affects aberration 
of the radio signal. At the same time, the Einstein de- 
lay is not affected by the finite emission height as Ae is 
simply the difference between the proper time and the 
coordinate time of pulsar, which is independent of the 
geometry and physics of emission. 

When calculating the timing contributions we will ne- 
glect the effect of gravitational lensing by the companion 
which mainly affects the Shapiro delay and may be rele- 
vant only for highly inclined binary systems |l8l Rlj, 123, 
l2lj . We will also neglect the retardation effect [221 123| 
which can be easily accounted for afterwards. 



A. Romer Delay 

We begin with the Romer delay Ar, the variation of 
the photon travel time as the pulsar moves through its 



A 



/,> - — no • R e 

c 



(15) 



where R e is the position of the emission region relative 
to the binary barycenter. The conventional DD model 
assumes R e = R p , while in our case R e = R p + hr, 
which gives rise to a correction to the Romer delay Ar^- 
We use equation ((THf to calculate Ah,/,, in which fo and 
A s r are constant and would not contribute to the timing 
signal, while A so r is small compared to A Q r by ~ e s (see 
below). Thus, we only retain A Q r and find 



A R<h 



(f3 OJ _ ■ m ) h 
3 cos 9 c 



ho 



3c cos 9 

HPo • y) 



[± (/3 D ■ x) sin $o sin a 



cos a — cos ( cos p 



sinC 



(16) 



where for any vector A, we denote its component per- 
pendicular to n as A_l = A — (A ■ n )n . Note that the 
first term in the right-hand side of this expression has 
a sign ambiguity, which reflects the term's opposite tim- 
ing contribution to the leading and trailing components 
of the pulse profile (corresponding to the two crossings 
of the same emission cone by vector d). This ambigu- 
ity has the effect of periodically changing the separation 
between the two components of the pulse profile on an 
orbital timescale, which makes this timing contribution 
similar to the so-called latitudinal time delay caused by 
orbital aberration 0, 0] • The second term in (|l(j[) ex- 
hibits no such sign ambiguity, meaning that it is the same 
for both pulse components, and manifests itself as an 
overall homogeneous shift of the pulse profile, similar to 
other conventional time delays. 

Assuming that all the angles specifying pulsar spin- 
magnetic orientation in (|16(l are of order unity (typical 
for rapidly spinning millisecond pulsars) one finds that 
the typical amplitude of Arj l is 



h a 

th — — Po 

c 



0.3 [is h 7 (/3 o /10~ 3 ) , 



(17) 



which is within the timing accuracy of the most stable 
pulsars. For slower pulsars which typically have small 
9 ~ p and the expression in brackets in (|16fl would 
be of order 9 <C 1, so that Ar^ ~ 9th, making this 
contribution difficult to measure. 

One may ask whether it is legitimate to discard A so r 
contribution in l|16|l when 9 -C 1. We believe it is, since 
the inclusion of |A so r| would add terms of relative sig- 
nificance ~ e s . At the same time, angles p and 9 are 
ultimately related to the position of the last open flux 
surface which touches the light cylinder at r^c — c/Q p . 
At the distance h this flux surface is characterized by an- 
gle 9 (h/rLc) 1 ^ 2 — £~J 2 ■ Thus, |A so r| contributions 
are smaller than Ar /, ~ 9th by 
neglected. 



e X J 2 <C 1 and can be 
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B. Shapiro Delay 

The Shapiro delay is caused by the curvature of space- 
time near the companion, which is only substantial when 
the pulse passes near the companion. It is accurately 
measurable only in nearly edge-on systems. The conven- 
tional formula for the Shapiro delay is 



A s ,o 



R 



log[n -(Rp-R c )-|Rp-R c |]+const, (18) 



where R g = 2GM c /c 2 is the Schwarzschild radius of the 
companion. 

A displacement of the emission region from the center 
of pulsar changes the propagation path of the radio beam 
and affects the Shapiro delay. The modified formula for 
the delay can be obtained by changing R p in equation 
(|18|l to R e = R p + hi and reads 



R u 



log [n • (R P + hr - R c ) 



|Rp + ftf-R c |] 
R 

— ^ log [(Rp - R c ) • (fi - n R ) 
hr ■ {h - n R )] , 



(19) 



where n^ = (R/, — R c )/|Rf, — R c |, and in the second line 
we have expanded the expression inside the logarithm in 
terms of h/\Rf> — R c |, while also setting f = fo. This can 
be rewritten as 



6b/\Rb — R c |. Thus, for millisecond pulsars with p, 8 ~ 1 
Aj,/, goes up to ~ (h/b)(R g /c) at conjunction. This 
can be substantial in a DNS binary: for a = 10 11 cm, 
h ~ 10 2 km, and i = 89° one has b m 17, 000 km so that 
h/b ~ 10 -2 and at conjunction Asm c & n reach ~ 0.05 fis. 

We also remark that the sign ambiguity of the first 
term in equation ill' 21) gives rise to somewhat different 
height-related contributions to the Shapiro delay for the 
leading and trailing edges of the emission cone, if their 
pulse arrival times are measured separately. 



C. Aberration 

Another time delay arises when the aberration from 
the velocity of the emitting region changes the spin phase 
of the pulsar at which its radio emission reaches the ob- 
server on Earth. In a conventional approach, aberration 
is due to /3 , while in our case it is due to (3 + /3 S . 

As demonstrated in [24[, DD, and [2lJ even in its con- 
ventional form (i.e. when (3 S = 0) aberration leads not 
only to the overall homogeneous time delay of the profile 
("longitudinal" aberration delay) but also to the distor- 
tion of the observed pulse shape. The latter can be inter- 
preted as the so-called "latitudinal" timing delay which 
has opposite signs for the different components of the 
pulse profile. The same result occurs in our case as well 
and we now consider these two effects separately. 



As 

Asj, = - 



A s ,o + Asm, 

R g hr ■ (n - n R ) 



(R 



p 



R c ) • (n - n R ) ' 



(20) 
(21) 



Far from conjunction the height related correction is 
very small: As.h ~ (h/ 'a)(R g / 'c) < 0.1 ns since h/a ~ 
10~ 4 and Rg/c ~ 10 /is even in DNS binaries. However, 
around conjunction (Rb — R c ) • (n — h R ) ~ 6 2 /|Rb — 
R c | -C a, where b ~ acosi is the minimum projected 
separation between the pulsar and its companion in the 
plane of the sky. At the same time, 

sin <&o sin a 

r -(n -n R ) = ± — — (n - n R ) x (22) 

6 cos v 

cos a — cos C cos p „ 

-{n -n R ) y 



3 cos 9 sin ( 
cos p + [3 cos 2 9 + 1] 5 
3 cos 9 



(n - n R ) z , 



and at conjunction (n — n R ) z ~ 6 2 /|Rb — R c |, while 
the first two terms in the right-hand side are of order 



1. Longitudinal shift 

The change A$ of the pulsar spin phase that arises 
as a result of aberration can be expressed using simple 
geometric arguments as: 



Ad ■ (s p x ftp) 
1 8 P x n | 2 



(23) 



Then, substituting in (|A12li - (IA14|) . the longitudinal tim- 
ing delay is found to be 



where 



a long _ A<I> _ .long . a long 
A - - ^A,0 + A,h ) 



^long _ Po ■ (Sp X fig) 

A0 ^ n p \s p X n | 2 



(24) 



(25) 



is the conventional longitudinal aberration delay [24J and 
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a long _ g s 
A,h ~ r, ■ 2 /■ 

i lp sin (, 

~w- -x Po± ■ m o cosC(3 cos 2 8 -\- l) 1 / 2 . „ . „ . r ,_ „ . „ .,"1 ,„„ N 

x J sin 2 C(/3 • *)(?„ • x + ° ° t — (ip • Po) - no • & n • r - cosC § P ■ f ) \ ■ 26 

^ 3 cost^ 3cos# J 



Here we have retained only terms variable on the or- 
bital time scale (proportional to O , i.e. orbital and spin- 
orbital contributions) . 

Orbital aberration alone gives A^ s ~ 1 — 3.3/xs. 
The height-related correction is smaller than by 
e s so that A^ 11 * ~ /3 D /ic _1 = th, where th is a fiducial 
timescale, see equation l)17|l . We discuss the measurabil- 
ity of the longitudinal aberration delay in sec. IIVI 

Note that unlike A^" S which is the same for all pulse 
components, the longitudinal height-related contribution 
AT^ contains terms with sign ambiguity (terms propor- 
tional to mo^ or ro, x ) 5 meaning that in the case of finite 
h even longitudinal aberration leads to pulse distortion. 



is the conventional latitudinal shift due to orbital velocity 
[24l | and A(am = (A so d • s p )/ sin£ is the contribution 
due to the finite height effects; A^A,h 1S different for the 
leading and trailing pulse components because of the sign 
ambiguity in mo >x and tq iX . 

Because of the pulsar beam shape, any shift in lati- 
tude results in the time delay of different components of 
the pulse profile (see Fig. 4 of Rafikov and Lai [2l| for 
an illustration). Leading and trailing components corre- 
sponding to the same emission cone are delayed by dif- 
ferent amounts 1211: 



2. Latitudinal shift 

Aberration also changes the latitude £ of vector d with 
respect to the pulsar spin axis. This latitudinal shift is 



Alat 
A 



±- 



' n p tan xo 



a lat 



a lat 



(29) 



ACa = , Ad '! P , = AC4.0 + A( A . h , (27) 
s p x n 



where 



A( A ,0 



Po±_ ' 

\s p x n | 



(28) 



where Ajf*^ is the latitudinal delay due to orbital aberra- 
tion (see eqs. [22]- [23] in |2l|) and the contribution due 
to the finite height effects is 



A 



lat 
A,h 



flp tan xo 



(f • x) [4cosC(n • p ) - (s p ■ p )] + 2 



{flo ■ {s p x f )) cosC(3cos 2 6» + l) 1 / 2 



(Po ■ X) 



tan £ 3 cos 9 

_ (Po± ■ m ) [(m ■ fig) - cosC(s p • m )] 
3 cos 9 tan £(fio ■ (s p x mo)) 



(30) 



Here yq is the position angle of linear polarization given 
by m 



cos a sin C — cos "l>o sm a cos £ 

The magnitude of the latitudinal shift is in general (as- 
suming tanxo ~ 1) of the same order as the longitudinal 
one, i.e. one should expect Ajf^ ~ th- However, in some 
cases, when the observer's line of sight cuts the emis- 
sion cone very close to its edge, tan \o may become very 
small which increases A^ ~ th/ tan^o, facilitating the 



detection of this delay [2l| . 

The height-related contribution A*f\ contains both 
terms which are the same for all pulse profile compo- 
nents and terms which have sign ambiguity. This is dif- 
ferent from the orbital latitudinal delay, which has iden- 
tical magnitude but opposite signs for different crossings 
of the emission cone. Based on the results of sec. IIII C II 
and IIII C 21 we conclude that height-related corrections 
change the simple separation of the effects of longitudi- 
nal and latitudinal delays, i.e. a homogeneous shift of 
the pulse profile in the former case and pulse distortion 



8 



in the latter. Both A^ s and AJI'q show a homogeneous 
shift as well as the distortion of the pulse profile. 



IV. MEASURABILITY 

Having calculated corrections to different time delays 
arising because of the finite height of emission region 
we can now ask whether these effects are measurable. 
The magnitude of these corrections depends on the type 
and orbital parameters of the system in question. In 
the double neutron star binary J0737-3039 @] one finds 
i/j ~ 0.33ft-7 /is, while the correction to the Shapiro delay 
is about O.O6/17 ns, (neglecting factors which depend on 
angles a, £, 77, etc.). At the same time, the current TOA 
uncertainty in this system is about 18 /js for a 30-s inte- 
gration at 820 MHz [2q[. which makes height-dependent 
corrections unobservable on short timescales in this sys- 
tem, despite its large value of f} = 10~ 3 . 

Much better timing accuracies have been achieved in 
some NS-WD binaries, however small companion masses 
make [3 significantly lower in these systems than in the 
DNS binaries. In particular, using the data from jj] we 
find t h w 21h 7 ns for J0437-4715 which has RMS tim- 
ing residual ctrms — 200 ns in 1-hour integration. For 
J1713+0747 we find t h « 12h 7 ns (a RM s = 125 ns in 
1-hour measurement), while for J1909-3744 th ~ 30hj 
ns (<7rms — 150 ns in 1-hour integration). Thus, fi- 
nite emission height effects are not currently observable 
in these systems. However, even more precise timing of 
these WD-NS binaries should allow future measurement 
of these effects by using many independent TOA mea- 
surements taken on longer time intervals. The ongoing 
projects such as PPTA @, and the advent of SKA 
all of which dictate understanding the systematic timing 
errors at the level of 1 ns, will enable the measurement 
of the finite height effects both by lowering the timing 
errors in already known systems and by discovering new 
binary pulsars which may be better suited for this goal. 

Another aspect of measurability has to do with the pos- 
sible degeneracy between the finite height corrections and 
the already known delays, which could arise if the orbital 
dependence of these corrections mimics the orbital vari- 
ation of other delays. In particular, DD have previously 
demonstrated that the longitudinal aberration delay 1 



A^f = ASty) + BCty), 

flbCLp sin r) 



cos i fl p cVl - e 2 sin £ ' 

S (ip) = sin ip + e sin ui, 
C(ip) = cosip + ecosu. 



(32) 
(33) 

(34) 



(with fit, = 2ir/Pi,) can be completely absorbed into the 
Romer delay by slightly rescaling the observed semimajor 



axis 
e.g. 



and eccentricity of the system from their true values, 



-.obs 



„obs 



cA 



(35) 



where a p = aM c /(M p + M c ). This makes A^ n s com- 
pletely unobservable on an orbital timescale. 

In fact, DD's argument can be generalized for any tim- 
ing contribution that depends on the pulsar orbital phase 
as XC{4>) + YS(ip), where X and Y are constant coef- 
ficients. The part of such a delay that is the same for 
all pulse components is then degenerate with the Romer 
delay. By the same argument, the sign-ambiguous de- 
lay contributions (which are different for the leading and 
trailing pulse features and cause pulse profile distortions), 
are degenerate with the conventional latitudinal delay 
AjfV which has the orbital dependence analogous to 
(JsSJl, see Rafikov and Lai |2lj . 

Using equations l(TB|) . and i(3Tj|) we see that A^, 
Ajwf j an d ^ Ah are Ihiear functions of /3 , which in the 
coordinate system (x, y, z) is given by 



00 



= {[cosi cos r]C(tp) + sin^S^)] x 



cVT 

[cos i smr]C(ip) — cosr/S^jy 
s'miC(ip) z}. 



(36) 



Thus, these finite height contributions are of the form 
(13211 with their own constant coefficients of order th- 
Capitalizing on the result of DD we can immediately 
say that all these height-dependent contributions are de- 
generate with A/j o an d A^q and are unobservable on 

an orbital timescale. The existence of A^ and A^ 1 ^ 
changes A and B in (|33|l to A + Ah and B + Bh, where 
Ah/ A ~ Bh/B ~ e s and the dependence of Ah and Bh on 
C, 77, i can be explicitly written out using equations 11GL 
(|26|l . and l|3U|l . Based on this argument we can also claim 
that the height dependent delays affect the observed val- 
ues of orbital parameters by changing ea — * £a + £h, 
where Eh ~ cth/a p = (h/a p )f3 ~ 10~ 7 . 

The orbital dependence of the Shapiro delay is unique 
enough to preclude As,h from being absorbed into cither 
Romer or latitudinal aberration delays. Nevertheless, we 
can look for a "new" orientation of the binary system 
characterized by the amended line of sight vector n new = 
h + Sn, where Sn is small and n • Sn = 0. Substituting 
fio = h new — Sn into i|19|) and choosing 



Sn — 



\Rp — R c 



rrcu 



(37) 



we find that 



1 This expression can be obtained by plugging 1361 into 1251 . 



R 

A s = — -5 log [n neu) • (R p - R c ) - |R P - R c |] , (38) 

i.e. the expression for the Shapiro delay reduces to its 
standard form l|18|) • In general the small correction vec- 
tor <5n varies with the pulsar orbital phase. However, the 
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Ag ; o attains its highest value near conjunction where the 
relative contribution of the Aj^ to A 5 is also largest 
and where |R P — R c | is roughly constant. Based on this 
we conclude that Ag^ can be incorporated into A$fl by 
slight reorientation of the binary system from no to n n ew ■ 
This makes Asm unobservable on an orbital timescale 
due to its degeneracy with A$fi- 

Despite the fact that timing contributions due to fi- 
nite height cannot be separately measured on an orbital 
timescale, one can hope to detect them on a longer geode- 
tic precession timescale, analogous to the measurement 
of the conventional aberration delay A^ s [DP: l2ll |2?| . 
Geodetic precession changes a pulsar's spin orientation 
which, leads to periodic variations of angles C and 77, af- 
fecting both the conventional aberration coefficients A 
and B and the height-dependent timing contributions. 
According to equation 1)35(1 . this would lead to periodic 
variation of the inferred values of the binary orbital pa- 
rameters (such as a and e), and would allow one to mea- 
sure the combined effect of orbital aberration and finite 
height effects. However, separating these two contribu- 
tors to the secular timing effects from each other may not 
be easy, as they vary in a similar fashion. In addition, 
the contribution to e from the finite emission height is 
rather small compared to sa- Additional information on 
the role of finite height effects compared to the orbital 
aberration may be provided by the pulse profile distor- 
tions on an orbital timescale (contributed both by A^'q 
and finite emission height effects), and by variation of the 
distortion pattern on the geodetic precession timescale. 

Probably the easiest and most promising way of de- 
tecting the finite emission height effects is by using ob- 
servations taken at different frequencies 2 . This method is 
based on the fact that the radio emission of different fre- 
quencies is produced at different heights in pulsar magne- 
tospheres [2, 113 1 a manifestation of the so called radius- 
to-frequency mapping (RFM) that will be discussed in 
more detail in sec. IV Al In this case h = h(v) and by 
using timing measurements at different frequencies one 
would infer a different value of a obs and e obs at each fre- 
quency since the height-dependent Eh is a function of v as 
well. As the orbital aberration is not frequency depen- 
dent, all variations of a obs and e obs that are frequency 
dependent must be ascribed to the finite emission height 
effects. By doing observations at sufficiently widely sep- 
arated frequencies, one can achieve corresponding height 
difference Ah ~ h, in which case the relative offset be- 
tween the values of a obs and e obs at different frequencies 
would be of order Eh ~ 10~ 7 . This estimate suggests the 
level of accuracy to which a obs and e obs must be measured 
at a set of frequencies in order for the height-dependent 
effects to be detected. 



2 This idea was suggested to us by Marten van Kerkwijk 



V. DISCUSSION 

Our use of a specific geometric model of the radio emis- 
sion (similar to RVM) does not imply that our results 
would be vitiated if pulsar beam has a non-circular shape. 
The assumption of a specific beam shape was necessary 
in our case only to carry out an explicit calculation of 
various height related delays in terms of angles 77, C, 9, 
etc. Clearly, even if pulsar has a non-circular beam, a 
finite emission height would still lead to aberration due 
to the pulsar spin and to the displacement of the emis- 
sion region from the pulsar center. Of course, in this 
case we would not be able to predict exactly how the 
height-dependent timing signals scale with the pulsar or- 
bital phase, which would complicate the interpretation of 
the secular evolution of these timing contributions on the 
geodetic precession timescale in terms of the pulsar spin 
orientation. However, even in this case one should still 
expect to see the difference in the observed orbital pa- 
rameters of the binary measured at different frequencies, 
see the discussion in the end of sec. IIVI 

We considered only the emission produced by a cir- 
cular source (tied to a single flux surface) which has the 
form of two 5-like spikes. This can be trivially generalized 
to the more realistic situation in which a set of different 
flux surfaces characterized by different values of 9 and 
p generates emission of different intensities which is ob- 
served on Earth as a complex pulse profile. This is done 
by a directly summing such signals with the proper delay 
for each flux surface, reflecting not only the conventional 
time delays but also the height-dependent contributions. 
Since the latter are explicit functions of 9, p, etc. which 
are different for different flux surfaces, it is clear that the 
finite height effects would lead to the periodic distortions 
of the whole pulse profile (even within only the leading or 
trailing part of the pulse profile) on an orbital timescale, 
as delays of different flux surfaces would vary differently. 
The relative temporal distortion of a pulse feature with 
width A<I> due to the finite height effects should be at 
the level of ~ A$fl p th which can reach ~ 1CP 4 for broad 
(A$ ~ 0.1) features in the profiles of the fastest (P p « 2 
ms) millisecond pulsars. 

In this respect we would like to emphasize the impor- 
tance of timing of the different components of the pulse 
profile separately rather than just fitting single template 
to the whole profile as is usually done. This procedure 
has been previously proposed for detecting orbital aber- 
ration on an orbital timescale [DD: l2lll27j . and it is clear 
from the previous discussion that it may also be very 
useful for interpreting the timing effects due to the finite 
emission height. 

Our study was done in a monochromatic approxima- 
tion while real radio receivers always have finite band- 
width. We do not expect this to affect our conclusions 
in any significant way, as the broadband pulse profile 
can be obtained by simply convolving the monochromatic 
pulse shape of the pulsar with the window function of the 
receiver. As long as one uses non-overlapping (even if 
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broad) bands one would still see the difference in the ob- 
served orbital parameters of the binary in different bands, 
which would in principle be a clear signature of the finite 
emission height effects. 



A. Delays due to frequency dependent Emission 
Height 

Our calculations have explicitly assumed the height of 
the emission region h to be constant. In reality, observa- 
tions [Til IT^l llfil ] find evidence for a radius-to-frequency 
mapping in which radio photons of different frequencies 
are emitted at different heights, so that h = h(v). 

It is then easy to see that the orbital motion of pulsar 
should give rise to additional time delays proportional 
to dh/dv. Indeed, because of the orbital Doppler shift, 
a monochromatic receiver detects radio photons which 
were emitted at slightly different frequency {Av/v ~ 
/3 ~ 10~ 3 ) than in the stationary case. The h{v) relation 
then implies that the height of emission also changes peri- 
odically to reflect this Doppler shift as 3 Ah/h ~ Av jv ~ 
(3 (observations suggest that d\nh/d\wu is roughly con- 
stant). This clearly translates into another correction 
to the Romer delay which is about Ah/c ~ f3 D h/c, i.e. 
has the same magnitude as Ap> h and all other height- 
dependent delays calculated in this paper. 

Another effect arises because of the direct relationship 
between the distance from the pulsar center h and the 
magnetic latitude angle 9 for a given flux surface: 



h{9) 



(39) 



where r e is the equatorial radius of the flux surface, as- 
suming that the dipolar geometry holds even beyond the 
light cylinder. The same Doppler effect that causes varia- 
tions in h also causes the emission cone to widen at some 
orbital phases and narrow at others with the amplitude 



A9 = 



tan 9 Ah 



(40) 



which causes pulse profile distortions. Because of that 
the direction of emission also changes by 

a dd . „ 3sin#(f + cos#rh) . . 

A ed = -A9 = - (3cQ ^ + 1)3/2 j A0, (41) 

see equation (j2J. Similar to the latitudinal aberration 
delay 21] such a distortion can be interpreted in terms 
of the corresponding time delay ~ lAgdlfi" 1 which, ac- 
cording to |@nj and ijHJ, is of order Hy 1 ^ sin 2 9 = 
f3 h/(fl p r e ). But since pulsar radio emission is produced 
in the open field line region, r e > r^c = c/il p , this ad- 
ditional delay is < (3 h/c - the same order of magnitude 
as all other height-dependent delays. 



There are two reasons why we have not explicitly 
computed timing effects arising from RFM despite the 
fact that they are non-negligible compared to the con- 
tributions which we have already calculated. First, all 
of these effects are proportional to d In h/d In v while 
the height to frequency dependence is not very strong, 
\d\ah/ dhiv\ 0.2 — 0.3 Second, the introduction 

of these terms does not change the qualitative picture 
of the finite emission height effects on timing and our 
conclusions about their measurability, but it does make 
calculations much more cumbersome. Because of this, we 
decided to leave the careful calculation of these effects for 
a future work. 



VI. SUMMARY 

We have extended the standard calculation of timing 
delays in binary millisecond pulsar systems by taking into 
account the fact that the emission region does not spa- 
tially coincide with the physical center of the pulsar but 
is located at some height (tens to thousands of kilome- 
ters) above the neutron star surface. This very natural 
extension has two consequences. First, such a displace- 
ment changes the propagation (Romer and Shapiro) de- 
lays. Second, the large spin velocity of the high-altitude 
emission region affects the aberration of the radio sig- 
nal and contributes to the aberration delay. We have 
calculated the height-dependent corrections to these de- 
lays in the framework of a specific model of the pulsar 
emission geometry. The typical magnitude of corrections 
(h/c)(3 is of order several hundred nanoseconds in com- 
pact DNS systems (which is at the current precision level 
of a few most accurately timed millisecond pulsars in NS- 
WD binaries), and should be measurable in the future 
by projects like PPTA and facilities like SKA. Although 
these timing corrections are degenerate with other de- 
lays (Romer and aberration), they can be detected using 
multi-frequency observations and exploiting the RFM. 
Detection of these effects could yield a new measure of 
the emission height in pulsars, which would provide a 
consistency check on the estimates using other methods. 
Although our calculation assumes a specific model of the 
emission region geometry (that of a circular emission cone 
in a dipolar magnetic field), the main results of our anal- 
ysis are qualitatively independent of that model. 
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3 The linear character of the relation between Ah and Av is not 
affected by the variation of the gravitational redshift in the field 
of pulsar and the change of spin velocity of the emission region. 
The coefficient in the relation changes because of these additional 



effects but only at the level R g /h and j3 s correspondingly, where 
R g is the Schwarzschild radius of the pulsar. 
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APPENDIX A: 



CORRECTIONS TO VECTORS m 
AND f 



Here we compute changes in the vectors f and m that 
result from the orbital and spin motion of the emission 
region. The corrections we are looking for are small, typ- 
ical orders are e so <C e a <C e s <C 1 for the spin-orbital 
A so , orbital A G , spin A s , and static terms correspond- 
ingly (the relative size of each of these terms is given in 
equations (O-JHI)). Thus, we can employ perturbation 
theory to compute these corrections by consecutively go- 
ing to different orders in these small parameters. 

In the static limit we neglect any motion of the emis- 
sion region in which case equations @, and © 
supplemented with the condition |m| = 1 result in equa- 
tions 10 and (fTUf) . 



1. Contributions due to pulsar spin 

In the next order we retain only terms proportional 
to the spin velocity f3 s c of the emission region, where 
(3 S = e s (s p x fo). This motion aberrates the vector d, 
producing an additional contribution A s d to do = no so 
that instead of © we now have 



d + A s d = n - f3 s 



(Al) 



as the condition for the pulse to be emitted towards the 
observer. For any vector A, we denote its component 
perpendicular to n as = A — (A • n )n . 

From the definition we find that (assuming constant 
6, i.e. no RFM, see sec. IVAl : 



A.,d = 



3cos#A s r — A s m 
(Scos^ + l) 1 ^ ' 



(A2) 



while equations QJ, Q and condition |rh| = 1 lead to 
three more constraints: 



A s m • m = 0, 
A s m • s p = 0, 
A s m ■ f o = -A s rm . 



(A3) 



Combining equations (|A1|) . (|A2|) . and l|A4(l we obtain 
six equations for six unknowns (components of A s m and 
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A s r), which can be easily solved: 



A s m 
A.„r 



(/3 S± ■ m )(s p x rh ) 

n • (§ p x mo) 
(f} s± ■ m )(s p x m ) 
3 cos ^no • (§p x rho) 
(Scos^ + l) 1 / 2 



3 cos 9 



(A4) 



(A5) 



providing us with the sought corrections due to the spin 
motion of the emission region. Note that although e 2 s can 
easily exceed e we do not need to go to higher order in 
e s because the spin contributions A s do not vary with 
the orbital phase of the pulsar and thus are undetectable 
by timing. On the other hand, despite the fact that the 
constant A s m and A s r given by l|A4|l and i|A5|l do not 
affect timing directly, we need to compute them in or- 
der to calculate the spin-orbital corrections in the next 
section. 



Substituting d = do + A s d + A D d in the right-hand side 
of this equation (in doing that we ignore the spin-orbital 
term A so d as it would lead to higher order corrections) 
one obtains 



= -0„± (A7) 
= -s s [f3 o j_(n ■ (§p x r )) + (n • (3 )(s p x r )± 
+ 2n (/3 o _L • (§p x f )) + (§p x A c r)j_]. (A8) 



A d 
A so d 



Finally, instead of l|A4j) we have the following relations: 

(A m + A so m) • s p = 0, (A9) 
(A m + A so m) • (f + A s r) = 

-(rh + A s m)- (A r + A so r), (A10) 

(A m + A so m) • (rh + A s m) = 0, (All) 



2. Corrections due to orbital motion and spin-orbit 
coupling 

In general the aberration of vector d is caused by the 
combination of both orbital and spin motions of the emis- 
sion region. This gives rise to additional orbital contri- 
butions A Q (which exist even if pulsar is not spinning or 
if h = 0, an assumption used by DD), and spin-orbital 
contributions A so to vectors rh and f , related to corre- 
sponding A D d and A so d via equations analogous to i|A2() . 

Instead of equation (|A1|) we now have 

d = d + A s d + A d + A so d 



which provide us with enough information to solve for 
A Q r, A Q m, A so r, A so m by considering separately terms 
of order e and e so . 

After tedious but straightforward calculations we find 



A Q m 
A r 



Sp x m 



n ■ (§p x m 
§p x mp 
3cos0rio • (§p x m ) 
(Scos^+l) 1 / 2 
3cos# 



09 o± • m ) (A12) 
{(3 0± ■ m ) 
Po±. (A13) 



iio — Po + d(/3 D • d) — j3 s + d((3 s ■ d). (A6) for the orbital terms and 



A.,„d 



(3 cos 2 8 + l) 1 / 2 

f3 o± (h ■ (ip x r )) + 2h (/3 o± ■ (s p x f )) (s p x [3 0± )± 

O COS O 



+ (n • o )(sp x f )_L + 



(§p x (§p x rh ))± P j_ ■ m 



3 cos( 



n ■ (§p x rh ) 



(A14) 



r 



for contributions to d resulting from the spin-orbital cou- separately in the timing calculations anyway). These ex- 
pling (we do not separately display expressions for A so r pressions are then used in sec. IHII to compute the height- 
and A so m as they are very cumbersome and are not used related corrections to different time delays. 



